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The unexpectedly strong suppression of high pt heavy-quarks in heavy-ion collisions has given rise 
to the idea that partons propagating through a medium in addition to energy loss by induced radia- 
tion also undergo substantial energy loss due to elastic collisions. However, the precise magnitude of 
this elastic energy loss component is highly controversial. While it is for a parton inside a medium 
surprisingly difficult to define the difference between elastic and radiative processes rigorously, the 
main phenomenological difference is in the dependence of energy loss on in-medium pathlength: in 
a constant medium radiative energy loss is expected to grow quadratically with pathlength, elas- 
tic energy loss linearly. In this paper, we investigate a class of energy loss models with such a 
linear pathlength dependence and demonstrate that they are incompatible with measured data on 
hard hadronic back-to-back correlations where a substantial variation of pathlength is probed. This 
indicates that any elastic energy loss component has to be small. 

PACS numbers; 25.75.-q,25.75.Gz 



I. INTRODUCTION 

The energy loss of hard partons propagating through the 
soft medium created in heavy-ion coUisions has long been 
regarded as a promising tool to gain information on the 
medium density evolution [l|, [1,1^ 0, H, @|- Radiative 
energy loss, i.e. the idea that medium-induced radiation 
predominantly carries away energy from a hard parent 
parton has been rather successful in describing not only 
the hadronic nuclear suppression factor for central colli- 
sions but also the effects of changing medium geometry 
Q. Calculations within dynamical evolution models in 
various formalisms @, ^, [13] have improved on this result 
and show also agreement with measured hard back-to- 
back correlatio ns [ill [l^ . \\^ and the measured suppres- 
sion of protons [l4| . 

Neverteless, there are indications that a radiative energy 
loss picture fails to describe the suppression of heavy 
quarks as seen in the single-electron spectra measure- 
ments [1^ and ener gy l oss due to elastic collisions with 
the medium [l^, [TtI llSl. [Tol . [20t has been suggested as a 
possible solution to this problem. Such calculations in- 
dicate a large component of elastic energy loss also for 
light quarks and gluons. 

In vacuum, the distinction between an elastic and a radia- 
tive process is straightforward — if the number of asymp- 
totic out-states is larger than the number of asymptotic 
in-states, the process is radiative. However, for a parton 
propagating through a medium, no asymptotic out-states 
can be defined. In particular, radiated quanta do not 
need to be on-shell as long as they re-interact with the 
medium within a sufficiently short amount of timescale as 
set by the uncertainty principle. It becomes thus to some 
degree a matter of convention if a particular process is 
seen as the (elastic) exchange of a virtual parton between 
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hard parton and medium or as the (inelastic) radiation 
of a virtual parton from the parent where the radiated 
parton is subsequently absorbed by the medium. 
However, while there is no sharp conceptual distinction, 
there are two crucially different regimes: If a parton ra- 
diated from the parent is highly spacelike, the formation 
time for this process is very short. On the other hand, 
the decoherence time for near on-shell radiation is very 
long, giving rise to interference effects and LPM suppres- 
sion. Phcnomcnologically, the first regime leads (in a 
medium with constant density) to a linear dependence of 
energy loss on pathlength whereas in the latter regime a 
quadratic dependence appears. 

It is the purpose of this paper to study observable conse- 
quences of this difference in pathlength dependence. To 
this end, we first construct a class of models in which 
the energy loss has parametrically a linear pathlength 
dependence. Based on the requirement that the nuclear 
suppression factor Raa for central collisions should be 
described, we point out limits for the parameter space 
of this model class. Finally, we demonstrate with three 
different scenarios the consequences of linear pathlength 
dependence for observables which explicitly probe path- 
length dependence such as Raa vs. reaction plane or 
hard back-to-back dihadron correlations and discuss lim- 
its for the relative magnitude of an elastic energy loss 
component. 



II. MODELLING ENERGY LOSS 

Key quantity for the calculation of energy loss in a dy- 
namically evolving medium is the probability distribution 
P{AE; E)path for a parton with initial energy E to lose 
the energy AE for any given path through the medium 
(in general, this includes a discrete part accounting for 
the possibility that the parton escapes without energy 
loss) . Empirically, one finds that a strong dependence of 
P{AE; E)path on the initial energy E does not seem to 
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be favoured by the data [2l| , thus we will in the following 
approximate P{AE;E)path « P{^E)path- 
From the energy loss distribution given a single path, 
we can define the overlap-geometry averaged energy loss 
probability distribution (P(A£'))t^^ as 

2^ /'27T I'OC I'OC 

{P{AE))t^^^— / d(f> / dxo / dyoP{xo,yo)P{AE)path. 

^TT Jo J-oc J~oo 

(1) 

Here, is the angle of the outgoing parton with respect to 
the reaction plane and P(a;o, yo) is the probability density 
for finding a hard vertex at the transverse position ro = 
(xo,yo) given the collision impact parameter b. This 
quantity is given by the product of the nuclear profile 
functions as 



P{xo,yo) = 



_ TA(ro + b/2)r4(ro-b/2) 



TAA(b) 



(2) 



where the thickness function is given in terms of 
Woods-Saxon the nuclear density /9A(r, z) as TaIi") = 
/dzp(^(r,2;). 

If we are interested in more differential observables, the 
relevant averaging procedure is changed, for example 
Raa as a function of the reaction plane is obtained from 



dxo / dyoP{xo,yo)PiAE)path 



(3) 



or the back-to-back per-trigger-yield can be computed us- 
ing {P{AE,E))Tr which is found by replacing P{xo,yo) 
in Eq. ([T]) by PTr{xo,yo^p), the conditional probabil- 
ity density to find a vertex in the transverse plane from 
which a near side hadron in the relevant trigger range 
was produced. The latter quantity is best calculated in 
a Monte-Carlo (MC) framework as it also has explicit 
momentum dependence, for details see [T2j . 
We calculate the momentum spectrum of hard partons 
in leading order perturbative QCD (LO pQCD) (explicit 
expressions are given in [121] and references therein) . For 
reasonably hard momenta (in practice 6 GeV or more) 
it can be assumed that hadronization takes place out- 
side the medium. Then the medium-modified observables 
arise from a convolution of the pQCD parton spectrum 
with the suitably averaged energy loss probability density 
and a hadronization function. For example, the medium- 
modified perturbative production of hadrons at angle 
can be computed from {P{AE))^ and the partonic cross 

section — =^^^^-7 using 



da 



AA~*h+X 



E 



,AA^f+X 

UjU 



^{P{AE))^^D}Zhi^,fi'p) 
(4) 

with D^'Zi^{z, fij^) the fragmentation function for parton 
/ with momentum fraction z at scale (23l . [2^ . From 
this we compute the nuclear modification function Raa 



vs. reaction plane as 



Raa{pt-,V, 



dN\J dPrdydcf) 
TAAih)daPP/dPTdyd<j)' 



(5) 



Thus, all the information about the medium and the en- 
ergy loss model is contained in P{AE; E)path- For ra- 
diative energy loss, we have in past works obtained this 
quantity by evaluating the two line integrals 



and 



lo 



(gL>(ro,(/-)= / dCm) 
Jo 



(6) 



(7) 



along the path of the parton through the medium and 
using the results of ^ to convert the result to a proba- 
bility distribution where we have assumed [25l . [2^ 



^(0 •2-e3/''(^)(coshp(^)-sinhp(0cosa). (8) 

Both the local energy density e(^) and the transverse 
flow field p(^) (with a the angle between flow vector and 
parton propagation) have to be inferred from a dynamical 
evolution model. For ease of comparison with published 
results, we use a 2-d hydrodynamical evolution model 
as in [13] and a 3-d hydrodynamical model [28] as in 

In a static medium, q{^) = const, and thus Eq. ([6]) ex- 
hibits the quadratic pathlength dependence for the en- 
ergy scale parameter Uc of radiative energy loss. 
We do not attempt to construct a first-principles model 
of elastic energy loss in the following. Such an attempt 
faces a number of difhculties, among them the running 
of the strong coupling as into a soft regime, the treat- 
ment of interactions of the high-pT parton in the hadronic 
phase of the medium, the precise nature of scattering 
partners in the QGP phase and q g conversion re- 
actions. Instead, we try to construct a class of mod- 
els which has the expected parametric dependence of 
energy loss probability distributions on the pathlength 
and infer then from the data limits for the parameters 
which characterize the model class. Since any observable 
considered here involves massive geometrical averaging 
(Eqs. ([I]),®) and since furthermore the nuclear suppres- 
sion factor is not very sensitive to the shape of even the 
averaged {P{AE))taa [IHj we may safely assume that 
details in the modelling of P{AE)path will not influence 
the outcome significantly. 

Note that in an ideal quark-gluon plasma (QGP), 
e'^/'*(^)(cosh p(^) — sinhp(^) cos a) is a measure of the en- 
tropy density and hence (up to a numerical constant) of 
the number density pM of scattering centers which can be 
probed by a parton on its trajectory ^. If the cross sec- 
tion for elastic energy loss is aei and the medium density 
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Pm, the number of scatterings dN on a path is given 
by dN = apMdS,. Exponentiating this expression, we ex- 
pect that the escape probabihty Pq without undergoing 
energy loss for a parton is parametrically given by 



exp 



-const. ■ aei / PAiiOd^ 



exp [—7 • K\ 



(9) 



where we have assumed that a^i is approximately inde- 
pendent of ^ and K is defined in analogy with Eq. ^ 
as 



K = J dCe^/*(C)(coshp(0 -sinhp(Ocosa) 



(10) 



taking into account the flow corrections to the probed 
density. Here 7 is a parameter with dimensions of a cross 
section measuring the interaction strength, and hence 
7g = 9/4:jq must hold to account for the different color 
factors of quarks and gluons. 

Note that in the case of radiative energy loss as calculated 
using the assumption of [J, [ll] the escape probability is 
determined by the condition of no radiation beyond the 
vacuum shower evolution rather than by a no-scattering 
condition. It can then be cast into the form ^] 



rad 



lim exp 



duj ^ 



(11) 



which requires the knowledge of the spectrum of medium- 



induced radiation 1 



. dljuj) 



as a function of radiated energy 



uj. This in turn depends on a decoherence condition for 
radiated quanta as well as LPM interference and cannot 
easily be cast into a simple form comparable to Eq. ([5]). 
A goo d discussion of the underlying physics can be found 
in^. 

If the parton does not escape without energy loss, it must 
undergo a shift in energy (there is also the possibility 
that a strong shift into a thermal regime occurs, which is 
equivalent to an absorption of the parton). It is reason- 
able to assume that the mean value of the shift in energy 
will grow linear in the number of scatterings N as 

with AEi the mean energy loss per scattering whereas 
the fluctuations around the mean will grow like VTV. As- 
suming a Gaussian distribution, this leads to the ansatz 



P{AE)patk = PoS{AE) + AAexp 



{AE ~ anf 



(12) 



where A/" is a normalization such that P{AE) = 1 
and (|12|) has to hold for quarks and gluons separately 
due to the different color factor, a is a parameter with 
the dimensions of a cross section times the energy shift 
per reaction. 

This class of energy loss models is characterized by three 
parameters: 



• a controls the mean shift in energy per expected 
scattering 

• P governs the strength of fluctuations around this 
mean shift. If f3 is small, the model will have a 
strong correlation between path (and hence initial 
vertex) and shift in energy, if the parameter is large, 
this correlation is lessened 

• 7 finally determines the magnitude of the escape 
probability. 

In a microscopical model, the actual distribution would 
presumably not be strictly Gaussian and the three pa- 
rameters would be correlated and calculable. Here, we 
will however pursue a different approach and see what 
can be inferred from the data. 



III. DETERMINING THE PARAMETER 
RANGE 

Our aim is to determine parameters of the elastic en- 
ergy loss model such that Raa for central collisions is de- 
scribed well and to use these parameter settings to com- 
pute observables where a different geometrical avergaing 
is performed. It does not seem straightforward to choose 
from the model space given by all possible {a, (3,^) all 
the solutions which are compatible with Raa- However, 
according to '^I'l, at least one parameter is fixed by the 
normalization of Raa- 

Let us consider 7 first: Clearly, the averaged discrete 
escape probability (Pq) may not be larger than Raa- On 
the other hand, if it is smaller than the error on Raa, 
this term is irrelevant. Via 



/'27r /"oo /"oo 

-z- I d(j) I dxo f dyoP(a;o,2/o)exp[-7K] < Raa (13) 

this translates to (evolution-model dependent) con- 
straints on 7. 

Let us now consider a scenario which is chosen such that 
(Po) is large and close to Raa- It follows that there 
must be a constraint on the shape of the continuous shift 
probability P[AE)path- The essential idea is apparent 
from the following: 

Averaging over both parton species and neglecting the 
fragmentation (which is a subleading correction on the 
shape of Raa m the high px region fT3|), Raa can be 
estimated as 



Raa[pt) ~ (Po)t..+ / (P(Ai?)T., 

JO+c 



dk: 



^(PT + AP) 



dNr, 



(14) 

where the lower integration boundary does not include 
the discrete contribution at and E^ax is the kinematic 
limit for the parton energy. If the spectrum is approxi- 
mated by a power law ~ then 
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AE [GeV] AE [GeV] 

FIG. 1: (Color online) Left panel: Continuous energy loss probability distribution for a quark propagating from the medium 
center outward in three different scenarios for elastic energy loss (see text). Right panel: Collision-geometry averaged energy 
loss distribution for quarks for three different scenarios. 



dN. 
dk 



dkq 



1/(1 + 

Pt 



In other words, Raa can be written as a constant term 
plus a term which increases with px and corresponds 
to the integral of the energy loss probability density, 
weighted by a steeply falling spectrum. From this term, 
sizeable contributions to Raa will come if {P{/S.E)taa) 
contains a lot of strength close to AE = (in this case 
the weight factor is close to unity) or if the spectrum is 
flat, i.e. n is small (as at the LHC [1^) or li pr ^ AE, 
i.e. at high pT- Thus, there is the generic expectation of 
a rise of Raa with pt [lOl ■ 

If 7 is already close to the allowed limit, a and /? must 
thus arrange in such a way that the shift term is small 
in order not to violate the limit. This means that the 
continuous energy loss distribution cannot contain much 
strength close to AE — 0. 

There are two possible ways this could happen (note 
that the continuous part must integrate to approximately 
1 — Raa)'- First, the distribution could be very flat and 
extend to large AE, thus the contribution in any given 
fixed interval in AE' would be small. Or second, the dis- 
tribution could be comparatively narrow but peaked at 
some large AE. 

In the case of radiative energy loss, the first scenario is 
realized, and for RHIC kinematics {P{AE))taa extends 
to O(IOO) GeV However, this cannot be so in the 

elastic case. Assume that along some fixed parton path 
elastic and radiative processes lead to the same mean 
energy loss. Then, a comparison of Eqs. (O and (fTO)) 
shows that parametrically the radiative energy loss will 
be smaller for all shorther paths but larger for all longer 
path. Thus, the dynamical range of expected energy loss 
generated in the model between short pathlength contri- 
butions close to the surface and long paths traversing the 
whole medium is vastly greater for radiative energy loss, 
making it very difficult to obtain a flat distribution in 





Scenario 1 


Scenario 2 Scenarios 


a [GeV"'] 


0.35 


0.35 0.04 


P 


7.0 


0.07 6.4 


7 [GeV-2] 


0.085 


0.085 0.12 



TABLE I: Parameters for the three different elastic energy 
loss scenarios described in the text 



elastic energy loss unless (3 is set to (unnaturally) large 
values. However, when this is done, the probability dis- 
tribution even for a single path is very wide, i.e. large 
fluctuations destroy the position-energy loss correlation 
and hence the tomographic information. 

We investigate two different scenarios with a large dis- 
crete escape probability (chosen such that it is compara- 
ble with the discrete quenching weights in the radiative 
energy loss scenario [22|). The parameters of each sce- 
nario can be found in tableU In the first one (Scenario 1) 
we assume a large value of (3, i.e. sizeable fluctuations of 
the energy loss given a path, in the second one (Scenario 
2) we choose small (3. The continuous energy loss proba- 
bility distribution using the medium evolution provided 
by the 3-d hydrodynamics for quarks for a single path 
from the medium center and and the geometry-averaged 
distribution {P{AE))taa for both scenarios are shown in 

Fig.m 

We also investigate a third scenario (Scenario 3) in which 
the discrete escape probability is adjusted to half the 
value of scenario 1 and 2. By the same argument seen 
above, the integral term in Eq. (jl4p must contribute 
more, hence the probability distribution needs strength 
close to AE = 0, leading to some pr dependent growth 
of Raa- This can only be achieved by choosing a com- 
paratively large (3. The resulting change in the shape of 
Raa due to the increased rise with pT as compared with 
the other scenarios disfavours even larger values of 7 (or 
a complete absence of escape without energy loss). 
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IV. Raa vs. reaction plane 

In Fig. [2] we show the resulting Raa for different an- 
gles with respect to the reaction plane for central Au- 
Au collisions at 200 AGeV and for non-central collisions 
at impact parameter & = 7.5 fm based on the medium 
description of the 3-d hydrodymnamics code [1^. For 
comparison we also include the radiative energy loss cal- 
culation. 

By construction, differences between the scenarios are 
not very pronounced in central collisions. However, the 
extrapolation to non-central collisions is rather different 
in all cases. First, all three elastic energy loss scenar- 
ios exhibit a pronounced rise with pT which is stronger 
than in the radiative case. As outlined above, this is con- 
nected with the strength of {P{AE))^ close to AE — 0. 
However, second and possibly more impartant, the split- 
ting between in-plane and out-of-plane emission in elastic 
energy loss is at most half of what is seen for radiative 
energy loss. 

This is not an unecpected feature: In a dynamic medium, 
the leading initial density dilution goes as l/r due to 
the longitudinal expansion (changing into with 
c(T) > 1 and growing later as pressure gradients drive ac- 
celerated expansion of the medium) . Inserting this result 
int Eq. pO| shows that elastic energy loss is dominated 
by early times when the medium is dense. At late times, 
K increases only logarithmically. Using typical evolution 
models, this implies that energy loss models with a linear 
pathlength distribution lead to sizeable losses only in the 
first 2-3 fm/c evolution time and become insensitive to 
late time behaviour. However, if the medium becomes ef- 
fectively transparent after 2-3 fm, a parton cannot probe 
density gradients much larger than this. This implies a 
loss of the sensitivity to in-plane vs. out-of-plane emis- 
sion for some partons, and even more dramatic, a loss 
of sensitivity to the difference between average near-side 
and away side pathlength. Consistent with the expec- 
tation, the smallest splitting is observed when P is set 
to large values and position-energy loss correlations are 
weakened in addition. 



V. BACK-TO-BACK DIHADRON 
CORRELATIONS 

Using the MC code of we compute the strength of 
hard hadronic back-to-back correlations for all three sce- 
narios of elastic e nerg y loss and compare with the STAR 
measurement [s^ |33| 

As can easily be seen, elastic energy loss gives a para- 
metrically wrong result — the calculation overshoots the 
data by a factor of more than two. This expectation was 
already formulated in [Tlj: In the radiative energy loss 
calculation, the difference between near side emission and 
away side emission arises from the drastic differences in 
average pathlength on near and away side p^ . However, 
since in all elastic scenarios only pathlengths of 2-3 fm 



can be probed before the medium becomes effectively di- 
lute, such mechanism cannot play a role. Thus, the main 
difference between near side and away side is that partons 
on the near side initially propagate into a zone of lower 
density and partons on the away side into higher den- 
sity. This accounts for the fact that there is additional 
suppression on the away side, however it is by far not 
as strong as in the radiative case. If the position-energy 
loss correlation is further weakened (as in Scenario 3) the 
disagreement with the data is more pronounced. 

VI. DISCUSSION 

It seems clear that a scenario in which energy loss is ex- 
clusively elastic is not very realistic — in reality, energy 
loss is presumably caused by different component pro- 
cesses which have different pathlength dependence. We 
may however estimate an upper bound for the strength 
of an elastic component with a simple ansatz 

PiAE)path = f ■ P{AE),i + (1 - f)P{AE)rad (15) 

where a fraction / of the energy loss probability distri- 
bution is caused by an elastic channel. 
The biggest uncertainty in the calculation of back-to- 
back correlations is given by the d-Au baseline for the 
correlation strength. If this uncertainty is combined with 
the uncertainty of the yield per trigger in Au-Au col- 
lisions and scenario 1 is chosen as being closest to the 
data, the 4-6 GeV momentum bin provides the strongest 
constraint and marginally allows an elastic energy loss 
fraction / < 0.32. However, this assumes that fragmen- 
tation is the only source of hadrons in this momentum 
bin which is not consistent with the fact that both radia- 
tive and coUisional energy loss individually underpredict 
Raa in this momentum window (cf Fig. [J). Thus, the 
true upper limit on / is even smaller. On the other hand, 
the most likely value of / in the analysis comes out rather 
smah as / = 0.086. 

In modelling the elastic energy loss, we have made several 
simplifying assumptions: First, we have assumed eikonal 
propagation of partons whereas in reality elastic energy 
loss leads to a deflection of partons. Furthermore, in 
writing down Eq. (jl2p we have assumed that aei doesn't 
strongly depend on the medium. However, if as{T) (and 
hence the interaction strength) grows very strongly for 
T Tc this could to some extend effectively result in a 
deviations from a strictly linear pathlength dependence 
(the same effect would however occur for radiative en- 
ergy loss). Furthermore, we have assumed a Gaussian 
distribution for the elastic energy loss probability. We 
have checked that the results do not crucially depend on 
this point. The reason is that the substantial geometrical 
averaging over many paths, cf. Eqs. (|1I3[) erases any de- 
tailed information on the shape of the distribution given 
a single path. In particular, we verified that replacing 
a Gaussian shape by a box shape with the same r.m.s 
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FIG. 2: (Color online) Left panel: Calculated nuclear suppression factor for radiative energy loss [g] and for the three scenarios 
of elastic energy loss described in the text for central collisions (impact parameter 6 — 2.4 fm) as compared to the PHENIX 
data [sJl Right panel; As on left panel, except for b — 7.5 fm. In-plane emission is indicated by solid lines, out-of-plane emission 
by dahsed lines. 




FIG. 3: (Color online) Away side yield per trigger for central 
200 AGeV Au-Au collisions for a trigger range between 8 and 
15 GeV as a function of momentum bin as compared with 
STAR data [3^ . [3^ . Shown is radiative energy loss in com- 
parison with the three scenarios for elastic energy loss (see 
text) based on a 2-d hydrodynamical model for the medium 
evolution [27l ]. 

width does not alter the resuhing Raa by more than 5% 
and yields the same back-to-back correlations within the 
statistical errors of the MC simulation. 
Finally, we have neglected any energy dependence of the 
energy loss probability distributions and hence finite en- 
ergy corrections. However, while these (and other) effects 
would clearly have an effect on the energy loss probabili- 
ties (to details of which the result is largely insensitive) it 



is less clear how they could possibly compensate for the 
generically different pathlength dependence of the two 
regimes. 

ft appears that the general features of this analysis are 
rather robust and there is a physics reason why this 
should be so: The difference in near side and away side 
suppression is consistent with a quadratic pathlength de- 
pendence but not with a linear dependence of energy loss. 
As apparent from the rather consistent overprediction of 
the away side yield even for large discrete escape proba- 
bility (the best possible scenario), this is a quite generic 
statement which holds for a whole class of elastic en- 
ergy loss models and does not depend on details of how 
the probability distribution appears microscopically. It 
is hard to imagine a mechanism which would change this 
picture without strongly changing the parametric depen- 
dence on pathlength. This places rather stringent limits 
on the relative magnitude of elastic contributions to light 
quark and gluon energy loss which may help to constrain 
microscopical models of energy loss. 
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